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The evolution of the curvature perturbation is highly non-trivial for curvaton models 
with self-interactions and is very sensitive to the parameter values. The final perturbation 
depends also on the curvaton decay rate F. As a consequence, non-gaussianities can be 
greatly different from the purely quadratic case, even if the deviation is very small. Here 
we consider a class of polynomial curvaton potentials and discuss the dynamical behavior 
of the curvature perturbation. We point out that, for example, it is possible that the non- 
gaussianity parameter /nl — while (?nl is non-zero. In the case of a curvaton with mass 
m ~ O(l) TeV we show that one cannot ignore non-quadratic terms in the potential, and 
that only a self-interaction of the type Vint = a 8 /A/ 4 is consistent with various theoretical 
and observational constraints. Moreover, the curvaton decay rate should then be in the range 
r = 1(T 15 - 1CT 17 GeV. 

§1. Introduction 

In the curvaton mechanism^ primordial perturbations originate from quantum 
fluctuations of a light scalar field which gives a negligible contribution to the to- 
tal energy density during inflation. This field is called the curvaton a. Inflation 
is driven by another scalar, the inflaton eft, whose potential energy dominates the 
universe. After the end of inflation, the inflaton decays into radiation. If the infla- 
tionary scale is low enough, H* -C 10~ 5 y / e^, the density fluctuations of the radiation 
component are much below the observed amplitude Sp/p ~ 1CP 5 and the fluid is for 
practical purposes homogeneous. While the dominant radiation energy scales away 
as p r ~ a -4 , the curvaton contribution to the total energy density may increase and 
the initially negligible curvaton perturbations get imprinted into the metric. The 
standard adiabatic hot big bang era is recovered when the curvaton eventually de- 
cays and thermalizes with the existing radiation. The mechanism can be seen as a 
conversion of initial isocurvature perturbations into adiabatic ones and, depending 
on the parameters of the model, is capable of generating all of the observed pri- 
mordial perturbation. The scenario sketched above represents the simplest possible 
realization of the curvaton mechanism, and a wide range of different variations of 
the idea have been studied in the literature. For example, the inflaton perturbations 
need not be negligible, 2 - there could be several curvatonsp* the curvaton decay can 
result into residual isocurvature perturbations®'^ and inflation could be driven by 
some other mechanism than slowly rolling scalarsP 

It is however well known that the predictions of the curvaton model are quite 
sensitive to the form of the curvaton potential.® - ^ 1 In particular, even small devi- 
ations from the extensively studied quadratic potential can have a significant effect, 
at least when considering non-Gaussian effects.®® One can also encounter strong 
scale-dependence of the non-gaussianity parameters.^ 1 When the initial curvaton 
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field value lies far in the non-quadratic part of the potential, the non-linear nature 
of the evolution equation will in general result in a very rich structure of phenomena 
in the parameter space, as has been discussed in detail mP^'^ 
The simplest non-quadratic curvaton potential is given by 

1 rr n+4 

V = -~vr?a 2 ' + A-r-p — , (1-1) 

2 M n ' v ; 

where n is an even integer to keep the potential bounded from below, and the in- 
teraction term is suppressed by a cut-off scale M. For non-renormalizable operators 
n > 0, we set the cut-off scale to be the Planck scale M = Mp = 1, and the coupling 
to unity, A = 1. For the renormalizable quartic case n = 0, the coupling A can be 
treated free parameter. 

The potential is reasonably well motivated by generic theoretical argu- 

ments. Indeed, the curvaton should have interactions of some kind as it eventually 
must decay and produce Standard Model fields. The curvaton needs to be weakly 
interacting to keep the field light during inflation. This however only implies that 
the effective curvaton potential should be sufficiently flat in the vicinity of the field 
expectation value during inflation but does not a priori require the interaction terms 
in (jl-ip to be negligible. Moreover, as typically the inflationary energy scale is rel- 
atively high, the field can be displaced far from the origin and therefore feels the 
presence of the higher order terms in the potential. The interactions could arise ei- 
ther as pure curvaton self-interactions involving the curvaton field a alone, or more 
generically as effective terms due to curvaton couplings to other (heavy) degrees of 
freedom that have been integrated out. An example of a possible physical setup 
which could lead to (jl-ip is given by flat directions of supersymmetric models that 
have been suggested as curvaton candidates.^ These would lead to a potential 
of the form (|l-ip with typically a relatively large power for the non-renormalizable 
operator. 

The amplitude and non-gaussianity of the perturbation depend on the curvaton 
decay time. Here we assume for simplicity a perturbative curvaton decay character- 
ized by some effective decay width F (for non-perturbative decay, see^'^). 

When the interaction term dominates in (jl-ip , the curvaton oscillations start in a 
non-quadratic potential and the curvaton energy density always decreases faster than 
for a quadratic case. For non-renormalizable interactions, the decrease is even faster 
than the red-shifting of the background radiation and the curvaton contribution to 
the total energy density is decreasing at the beginning of oscillations. Consequently, 
the amplification of the curvaton component is less efficient than for a quadratic 
model. For the same values for m and r, the curvaton typically ends up being more 
subdominant at the time of its decay than in the quadratic case. 

Despite the subdominance, the curvaton scenario can yield the correct ampli- 
tude of primordial perturbations as the relative curvaton perturbation bo* jo* pro- 
duced during inflation can be much larger than 10~ 5 . For a quadratic model, it is 
well known that the curvaton should make up at least few per cents of the total 
energy density at the time of its decay in order not to generate too large non- 
gaussianities.®'^ This bound does not directly apply to the non-quadratic model 
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(|l-ip since the dynamics is much more complicated. Although the subdominant cur- 
vaton scenario implies relatively large perturbation 5a* /a*, the higher order terms 
in the perturbative expansion of curvature perturbation can be accidentally sup- 
pressedP'^CS 

§2. The curvature perturbation 

Let us adopt the 5N formalisrrPSl 1 '^ and assume that the curvature perturbation 
£ arises solely from the inflation generated perturbation of a single curvaton field. 
Then 

C(i, x) = N'(t, Q5cr* (x) + ~iV"(t 5 Q 8a* (x) 2 + );N'"(t, Q5a, (x) 3 • • • . (2-1) 

2 b 

Here N(t,t*) is the number of e-foldings from an initial spatially flat hypersurface 
with fixed scale factor a(t*) to a final hypersurface with fixed energy density p(t), 
evaluated using the FRW background equations. The final time t is some arbitrary 
time after the curvaton decay. The prime denotes a derivative with respect to the 
initial curvaton value a* . Here we take to be some time during inflation soon after 
all the cosmologically relevant modes have exited the horizon and assume that the 
curvaton perturbations 8a* are Gaussian at this time. The expansion (|2-ip is then 
of the form 

C(t,x) = Cg(t,x) + ^/ NL Cg(t,x) 2 + ^ 5NL Cg(t,x) 3 + --- . (2-2) 
where C g (£,x) is a Gaussian field and the non-linearity parameters are given by 

25 N "' / ,N 

5NL = 54 N« • (2 ' 4) 

Here we neglect all the scale dependence of the non-linearity parameters''^ With 
this assumption, and neglecting higher order perturbative corrections, the constants 
/nl and #nl measure the amplitudes of the three- and four-point correlators of £, 
respectively. Observationally, they can be extracted from the CMB bi- and trispec- 
tra. 

We assume the curvaton obeys slow roll dynamics during inflation and introduce 
a parameter r* to measure its contribution to the total energy density p at t*: 



Pa 
P 



W « 1 • (2-5) 



Here the inflationary scale H* is a free parameter, up to certain model dependent 
consistency conditions. Assuming inflation is driven by a slowly rolling inflaton field, 
we need to require <C lCP 5 -^ in order to make the inflaton contribution to ( 
negligible. In this setup we also need to adjust the slow-roll parameter e = —H*/H? 
determined by the inflaton dynamics, to give the correct spectral indexpSJ n — 1 
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2e — 27] aa . The curvaton contribution, r\ ac = V"(a*)/3H 2 , is typically negligible 
because of the sub dominance of the curvaton. The curvaton mass is required to be 
small, but the same also holds for the inflaton mass. 

After the end of inflation, one assumes that the inflaton decays completely into 
radiation and the universe becomes radiation dominated. The decay constant r 
accounts for the coupling between the radiation and the curvaton component. The 
evolution of the coupled system is then given by 

a + (3H + r)d + m 2 a + X(n + 4)a n+3 = (2-6) 
p\ = -4Hp T + Fa 2 (2-7) 
3H 2 = p T + Pa . (2-8) 

The initial conditions are given by p T = 3H 2 and p a = V(a*) = r*/(l — r*)p r 
specified at time t* corresponding to the end of inflation. We also set a = 0. Given 
the parameters n, A and m, which determine the potential (jl-ip . and the two initial 
conditions r* and H*, one can calculate N in (|2-1|) from this set of equations. To 
find the curvature perturbation, we set 5a* = H*/(27r) and compute ( = N(a* + 
5a*) — N(a*). For a given set of parameters, one then adjusts the decay width r so 
that the observed amplitude^ 1 £ ~ is obtained. 

One may treat r as a free parameter since we have not specified the curvaton 
couplings to other matter, in particular to the Standard Model fields. However, since 
the primordial perturbation is mostly adiabatic, the curvaton should decay before 
dark matter decouples in order not to produce isocurvature modes. Assuming a 
freeze-out temperature T ~ ?71dm/20 for an LSP type dark matter model with the 
LSP mass todm ~ 0(100) GeV, this translates to a rough bound 

r > 10~ 17 GeV . (2-9) 

While this bound could be relaxed in non-minimal constructions, let us here adopt 
(|2-9p for definiteness. 



§3. Small deviations from the quadratic form 



Let us first assume that the deviation from the quadratic from is small and write 
the curvaton potential (|l-ip as^'^l 1 



V(a) 



1 2 2 , » 4 (°\ n 

-m a + Xm — 

2 \mJ 



(3-1) 



where A is some coupling constant. It is also useful to define a parameter s which 
represents the size of the non-quadratic term relative to the quadratic one: 

n-2 



2A 



rn 



(3-2) 



Thus the larger s is, the larger is the contribution from the non-quadratic term. 
The curvature fluctuation can then be written, up to the third order, as^S 

2 



C = 5N 



2 oi 
—r 

3 a. 



osc 5a* + 
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where a osc is the value of the curvaton at the onset of its oscillation; abova^] 

?>Po 



4/5rad + 3p 



(3-3) 



(3-4) 



decay 



Notice that o-' osc /a osc = l/o* for the case of the quadratic potential. With this 
expression, we can write down the non-linearity parameter /nl as 



NL 



where we have defined 



4r 



5 



(1 + 5) 



5 
3 



•5r 



r/2 



(3-5) 



(3-6) 



with 5 = for a purely quadratic potential. Also notice that, although the curvaton 
scenario generally generates large non-gaussianity with /nl ^ 0(1)) the non-linearity 
parameter /nl can be very small in the presence of the non-linear evolution of the 
curvaton field,©' 1 ® which can render the term 1 + S ~ . 

The non-linearity parameter ^nl associated with the trispectrum can be written 

as 



9NL 
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°osc°osc 



cr. 



/3 



+ 35 



- (1 + S) + I (1 - 95) + lOr + 3r 2 



(3-7) 



As one can easily see, even if the non-linear evolution of a cancels to give a very 
small /nl, such a cancellation does not necessarily occur for #nl- Examples® of 
the behavior of the non-linearity parameters for various types of self-ineraction are 
depicted in FigfTJ 



§4. General results 



Let us now relax the assumption that the deviations from the quadratic case are 
small. The numerical solutions of the equations of motion (|2-6[) exhibit complicated 
behaviour not qualitatively present in the simplified analytical approximation. To 
demonstrate this, in Fig. [2] we plot AN as a function of H(t), or the inverse of time, 
for fixed initial values of H*, r* (see (|2-5p ) and r for two different curvaton mass 
values. For different masses the moment of transition from the non-renormalizable 
part of the potential to the quadratic part of the potential is different, and this affects 
strongly the final value of AN which is dictated by the duration of oscillations in the 
non-quadratic regime. Deep in the quadratic regime the oscillations become faster 

*' Note that in the literature there is much variation in the definition of r, as is also in the 
present paper. 
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Fig. 1. (Left) Plot of /nl as a function of n for several values of s. (Right) Plot of <7nl as a function 
of /nl for several values of s. Notice that /nl and n have one-to-one correspondence. In both 
panels, r — 0.01. 



and faster so that AN's evolution is given by a scaling lawP^ However, before the 
quadratic term in the potential starts to dominate, AN shows complicated oscillatory 
behaviour that can be only tracked numerically. 




10 18 10 16 10 14 10 12 10 10 10 s 

H 

Fig. 2. AN as a function of H , (time evolves from right to left), for the masses m = 10 -11 (lower 
curve) and m = 5 ■ 10~ n (upper curve). In both examples H* — 10 -6 , n — 4, F = 10 -18 and 
r, = 10- 10 . 



From Fig. [2]it is clear that as the field value oscillates in time, so does AN. In the 
non-quadratic regime AN oscillates with a large amplitude. If the transition to the 
quadratic regime is slow compared to the oscillations in the non-quadratic regime, 
the transition averages over several oscillations. As a consequence, the final value 



The self-interacting curvaton 



7 



of AN will be a non-oscillatory function of the model parameters. However, if the 
oscillation frequency in the non-quadratic potential is slow, and the transition to the 
quadratic oscillations is rapid, then the phase of the non-quadratic oscillation affects 
the final value of AN. If the parameters happen to be such that the transition 
to the quadratic regime occurs at a maximum of the oscillation, a relatively high 
value of AN freezes out. Similarly, if the transition occurs at a minimum of the 
oscillation cycle, the final value of AN will be much smaller. If the parameters 
governing the moment of transition, such as the curvaton mass m, are changed 
continuously, then the phase of the non-quadratic oscillation during the transition 
also changes continuously. In the space of the parameters this results in an oscillatory 
pattern in AN. This behaviour can be understood by observing that the curvaton 
energy density at the beginning of its oscillations in the non-quadratic part of the 
potential can not be expressed in terms of an amplitude of the envelope alone but 
also depends on the phase of the oscillation, or equivalently on both the field a and 
its time derivative &, in a non-trivial way. In effect, these act as two independent 
dynamical degrees of freedom. If the transition from the interaction dominated part 
to the quadratic region takes place at this stage, the initial variation of the curvaton 
value cr* can therefore translate in a non-trivial fashion into the final value of the 
curvature perturbation. 

For a potential with V ~ <r n+4 no oscillatory solutions exist if n > 6. This 
means that in the non-quadratic regime the curvaton merely decays and hence no 
oscillations in AN occur, we have scanned the parameter space (to, n, T) 

to find the regions that yield £ ~ 10 -5 and acceptably small non-gaussianity while 
being consistent with the slow-roll assumption. For fixed parameter values, the result 
can be mapped out as a region in the space of the inflation scale and the initial 
relative curvaton energy fraction, r*. We also bound H* from above by H* < 10~ 5 , 
in order to prevent the excessive production of primordial tensor modes and to keep 
the inflaton perturbation negligible. 

The experimental limits for /nl axe given by the WMAP 5-year data, 27 ' —9 < 
/nl < HI; we also require that -3.5 x 10 5 < #nl < 8.2 x 10 5 as given in: 29 ' The 
schematic outcome of the scan is depicted in Fig. [3j 

The observational limits for /nl and #nl constrain the allowed area in the very 
subdominant regions of the parameter space, depicted in Fig. [3] by the line a. Other 
constraints shown in Fig. [3] arise from the internal consistency of the self-interacting 
curvaton scenario. The bound b is obtained because otherwise the initial perturba- 
tions would be too small to produce the observed amplitude. The bound c reflects 
the requirement that the curvaton should be massless, or V" < H^, which is neces- 
sary for the generation of curvaton perturbations during inflation. Because of the 
subdominance of the curvaton, the realistic bound should arguably be a few orders 
of magnitude tighter. However, a change of an order of magnitude moves the actual 
cut by a very small amount in the log- log plots. Finally, the bound d guarantees the 
absence of the isocurvature modes in dark matter perturbations and corresponds to 
the limit on the curvaton decay width given in (|2-9|) . Detailed figures that show the 
allowed regions of and r* for different parameter values can be found inJ^J 
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Fig. 3. A schematic illustration of the different cuts limiting the allowed area in the parameter space 
of the inflation scale H* and the initial curvaton energy density relative to the total density r* . 



§5. TeV mass curvaton 



It is of particular interest to consider the self-interacting curvaton with a mass 
m — 1 TeV. Examples of such a curvaton could be found among the MSSM flat 
directions®} or light moduli fields of string theories. Let us therefore fix m = 1 TeV 
and demonstrate first that a simple non-interacting or quadratic form does not give 
rise to a consistent curvaton model, as discussed inP^ 

For a quadratic curvaton potential one can write the perturbation as 



<7* 



10" 



(5-1) 



where H*/a* gives the initial perturbation amplitude in the curvaton, and r e g is the 
efficiency factor that can be approximated quite well by the energy fraction at the 
curvaton decay P* 

Pa 



?~eff ~ r dc 



Pr + Pa 



decay 



(5-2) 



Relating <7* and r* from ^m 2 a%/3Mp l H% ~ r*, and noting that r dec < 1, we find the 
constraint on the initial curvaton energy fraction 



1 m 2 
r* < -- 



(5-3) 



pi 
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In the free curvaton case r^ec also determines non-gaussianity through the simple 
relation^ 1 /nl = 5/4r,j ec . Very roughly, observationally |/nl| < 100, which implies 
the constraint 



10- 



n> —eW,- (54) 

The limits (|5-3p and (|5-4p are well known. However, there is more. Since the 
observed perturbations are adiabatic to great accuracy, the curvaton must decay 
before dark matter decouples. For each set of the initial conditions, (il*, r*), there is a 
relation between r^ec i n (15-2[) and the effective decay constant T given by the fact that 
decay time is defined as H = T. Here we assume implicitly a perturbative curvaton 
decay, but T could stand for any effective inverse decay time and thus the following 
discussion should hold, at least roughly, also for a non-perturbative curvaton decay 
as discussed irP^} (note however that non-perturbative curvaton decay could turn 
out to be a source of a considerable non-gaussianity^S} ) . 

The exact evolution of the energy densities is difficult to solve analytically. How- 
ever, we can approximate the curvaton evolution by dividing it up to three phases: 



1. 



3. 



H , the curvaton is effectively massless, so the field value 



When V" = m 2 < 
stays constant, a = 

When V" = m 2 > H 2 , the curvaton oscillates in the quadratic potential, and 
thus its energy density approximately scales as p a oc a -3 . 
The curvaton oscillates until H = T. whence it decays. 



Solving the Friedmann equation for the regime where m > H then yields 



a* 




Hit 
H 



1 + 



Hi 



rn 



\[H 



■rn 



O (rl) . 



Using the above result we can solve for r* to find 



in 



\JmT 



6(^1 C 



Hi Hi _ 12 ( Mm ) n 

my/mT V m ' 



(5-5) 



We need to check whether, given the constraints discussed above, the self-interactions 
can be neglected if m ~ 1 TeV. Thus, adopting the form of the potential given in 
(jl-ip . in order for the quadratic assumption to be consistent, we should require that 



n+4 



-m 2 a 2 > 
2 Mpj 



(5-6) 



throughout the evolution. Since the energy density of the quadratic field decreases 
monotonously, it is sufficient to apply this requirement only for the initial conditions. 
Solving for r* such that the magnitudes of the quadratic and non-quadratic terms 
are equal, we find the condition 
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We have plotted this condition for n = 4 in figure H] as the diagonal dotted line. 
To the right of it, the non-quadratic term dominates initially. As can be seen in 
figure HI there is practically no allowed region in the parameter space where the 
quadratic assumption would even approximately apply. For smaller values of n, the 
self-interaction becomes important even for much smaller values of and r*, and 
thus, there is no quadratic regime left in the parameter space. 




10 4 10 6 10 8 10 1(> 10 12 10 14 

Fig. 4. Parameter space of the quadratic curvaton. r* must be above the lower horizontal dashed 
line to produce £ ~ 10 -5 (equation H5-3| l1 and below the blueupper horizontal dashed line to 
produce small enough /nl (equation (|5-4[l . Furthermore, r is constrained from above, and thus 
only the parameter space to the right of the black solid line is allowed (equation l|5-5|l V The 
green dotted line illustrates the equality of the mass term and a possible self-interaction term 
in the potential (equation (15-7[) ) for n = 4. For smaller values of n the line moves further to the 
left. To the right of the dotted line the self-interaction dominates, and thus practically in all of 
the allowed parameter space the self- interaction must be taken into account. 



We thus may conclude that even if the curvaton self-interactions were very weak, 
a purely quadratic potential would not be a consistent approximation for a mass 
m — 1 TeV; instead, the effects of the self-interactions need to be taken into account. 
These change the dynamics of the curvaton in a significant way. Moreover, as dis- 
cussed inj^ a scan of the parameter space reveals that only n = 4 potential with 
V ~ a s has any allowed parameter space. In addition, in order to obtain a correct 
perturbation amplitude, the decay width r should be in the range 
GeV. For most particle physics models, this would be a rather small decay width. 
We estimate^ 1 roughly that in the MSSM, where the non-zero curvaton background 
provides masses to other particles and hence gives rise to a kinematical blocking,^ 
one could obtain widths of the order r ~ 10~ 12 . However, a detailed and more 
proper calculation is required to settle the issue. 
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§6. Discussion 

It may appear surprising that even very small deviations from the quadratic 
form of the curvaton potential can affect the curvature perturbation in a significant 
way. However, one should bear in mind that the small curvature perturbation is 
really the difference of two large numbers. The number of e-folds generated during 
curvaton oscillations is typically N ~ 0(10), whereas the difference that gives rise 
to the non-gaussianity is AN < 10~ 8 . Since self- interactions imply non-linearities 
in the evolution of the curvaton field and in the number of e-folds N, one can un- 
derstand that even small changes can have profound effects in the difference AN. 
In particular, as discussed here, the non-gaussianities turn out to be quite different 
as compared with the simplest quadratic model. There the magnitude of /nl in. the 
limit r^ec <C 1 is determined by the curvaton energy density at the time of its decay, 
/nl ~ 1/^dcc- However, with self-interactions the prediction for /nl can significantly 
deviate from this simple estimate. Even if r^ ec <C 1, there exists regions in the pa- 
rameter space with |/nl| < This is because the value of /nl oscillates and 
changes its sign. Nevertheless, (7nl can then be very large and one has a rather non- 
trivial non-Gaussian statistics characterized by a large trispectrum and a vanishing 
bispectrum. Such a situation, discussed already in,« appears to be rather generic in 
self-interacting curvaton models, and is possible for a wide, albeit restricted, range 
of model parameters. Large non-gaussianities can be generated even if the curvaton 
dominates the energy density at the time of its decay. In general, in the presence of 
self-interactions the relative signs of /nl and <?nl and the functional relation between 
them are typically modified from the quadratic case. Thus the non-linearity param- 
eters taken together, in possible conjunction of other cosmological observables such 
as tensor perturbations, may offer the best prospects for constraining the physical 
properties of the curvaton. 

A TeV mass curvaton is a rather special case. An important constraint, valid also 
for higher mass curvatons, is that it has to decay before the CDM freeze-out. This, 
together with observational constraints, fixes the range of the initial conditions for 
the curvaton field which turn out to be such that the quadratic term in the curvaton 
potential cannot dominate over possible higher-order terms for the whole dynamical 
range. One finds 17 ' that the only viable curvaton potential that satisfies all the 
constraints is V = m?a 2 /2 + cr 8 /M 4 . Moreover, the curvaton decay rate should 
be in the range T = 10 -15 — 10 -17 GeV. Note that in the case where the curvaton 
energy density is subdominant at the time of decay, the curvaton does not necessarily 
have to decay before baryogenesis, which can be a process that takes place among 
the inflaton decay products. However, the decay should be able to produce thermal 
CDM particles so that the CDM perturbation is adiabatic. 

Note also that what really matters is the equation of state, not the time of decay. 
Thus if the curvaton decays too early, the perturbations might still generated if the 
decay products have the equation of state of matter. An example of this could be 
the MSSM flat direction fragmenting into Q-balls, which would then slowly decay. 
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